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Abstract 

Recently, it was predicted theoretically and verified experimentally that a pair of delayed and 
cross-polarized short laser pulses can create molecular ensembles with a well defined sense of 
rotation (clockwise or counterclockwise). Here we provide a comparative study of the classical 
and quantum aspects of the underlying mechanism for linear molecules and for symmetric tops, 
like benzene molecules, that were used for the first experimental demonstration of the effect. Very 
good quantitative agreement is found between the classical description of the process and the 
rigorous quantum mechanical analysis at the relevant experimental conditions. Both approaches 
predict the same optimal values for the delay between pulses and the angle between them, and 
deliver the same magnitude of the induced oriented angular momentum of the molecular ensemble. 
As expected, quantum and classical analysis substantially deviate when the delay between pulses is 
comparable with the period of quantum rotational revivals. However, time-averaged characteristics 
of the excited molecular ensemble are equally good described by the these two approaches. This 
is illustrated by calculating the anisotropic time-averaged angular distribution of the double-pulse 
excited molecules, which reflects persistent confinement of the molecular axes to the rotation plane 
defined by two polarization vectors of the pulses. 
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I. INTRODUCTION 



Laser control of molecular rotation, alignment and orientation has received significant 



attention in recent years (for a review, see e. g. pj, 121]). Interest in the field has increased, 
mainly due to the improved capabilities to manipulate characteristics of the laser pulses 
(such as time duration and temporal shape), which in turn leads to potential applications 
offered by controlling the angular distribution of molecules. Since the typical rotational 
motion is 'slow' (~10ps) with respect to the typical short pulse (~50fs), effective rotational 
control and manipulation are in reach. During the last decade, temporal rotational dynamics 



of pulse-excited molecules was studied 



IJ] , and multiple pulse sequences giving rise to the 
enhanced alignment were suggested Q-?] and realized jsH]. Further manipulations such 
as optical molecular centrifuge and alignment-dependent strong field ionization of molecules 
were demonstrated jl^ . Selective rotational excitation in bimolecular mixtures was 



suggested and demonstrated in the mixtures of molecular isotopes [14| and molecular spin 



isomers 



ll6|. 



Recently, several groups suggested a method for exciting field-free unidirectional molec- 
ular rotation, in which the angular momentum along the laser beam propagation direction 
is provided by two properly delayed ultrashort laser pulses that are linearly polarized at 45 
degrees to each other 17142 1|. The first experimental demonstration of the effect was done in 



2l| for benzene molecules. The mechanism behind this doub 



and is most easily explained in the case of linear molecules 17H19| . The first ultrashort laser 



e-pulse scheme is rather clear. 



pulse, linearly polarized along the z-axis, impulsively induces coherent molecular rotation 
that continues after the end of the pulse. The molecules rotate under field-free conditions 
until they reach an aligned state, in which the molecular axis with the highest polarizabil- 
ity is confined in a narrow cone around the polarization direction of the first pulse. The 
second short laser pulse is applied at the moment of the best alignment, and at an angle 
with respect to the first pulse. As a result, the aligned molecular ensemble experiences 
a torque causing molecular rotation in the plane defined by the two polarization vectors. 
The torque acting on a linear molecule is maximal when the laser pulse is polarized at 45 
degrees with respect to the molecular axis of the highest polarizability. This defines the 
optimal angle between the laser pulses. The direction of the excited rotation (clockwise or 
counter-clockwise) is determined by the sign of the relative angle (±45 degrees) between the 
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FIG. 1: (Color online) "Molecular propeller" (figure taken from reference 



first and the second pulse in the polarization plane. This double pulse scheme (see Fig. [T]) 



was termed "molecular propeller" |17H19j. as it resembles the action (side kick) needed to 



ignite a rotation of a plane propeller. The heuristic arguments presented above are, in fact, 
essentially classical, although it is clear that a detailed quantum analysis may be needed to 
consider the long-time evolution of the excited molecules, or the operation of this scheme 
at relatively weak pulses when only low- lying rotational states are involved. That is why 



all the previous treatments of the problem 



17 



2lj presented full-scale quantum-mechanical 



studies, combining both analytical and numerical methods. 

In the present paper, we provide a detailed comparison between the classical and the 
quantum-mechanical analysis of the double pulse scheme for exciting unidirectional molecu- 
lar rotation, both for linear molecules and for benzene (that was the object of the experiment 



m 



2l|). We will demonstrate the predictive power of the classical treatment that correctly 



prescribes the optimal delay between the pulses, and the angle between them, and will 
find the advantages and limitations of such an approach in quantitative description of the 
effect at experimental conditions. We believe that by comparing the complementary clas- 
sical and fully quantum mechanical views on the same phenomenon, the reader will get a 
comprehensive idea about the underlying physics. 



Recent paper 



19| studied quantum-mechanically the "molecular propeller" scheme in 



the case of linear molecules. Dependence of the induced angular momentum on the pulses' 
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intensity, and the delay between them, was analyzed. Moreover, the angular distribution of 
linear molecules subject to unidirectional rotation was considered, and it was shown to be 
confined in the plane defined by the two polarization vectors of the pulses. This anisotropic 
angular distribution was characterized by the observable (cos^y?), which was referred to as 
the azimuthal factor. 

In the present paper, we describe the same scheme using a classical description of the 
double-pulse molecular excitation, followed by a Monte Carlo averaging over the thermal 
ensemble. We compare the time dependence of the azimuthal factor according to the classi- 
cal analysis, to the quantum mechanical one, and show that both results coincide for initial 
times. In addition, we introduce a long-time-averaged angular distribution function, which 
conveniently depicts the molecular orientation following the excitation. This distribution 
function seems to behave similarly in the classical and quantum cases. Finally, we analyze 
classically the excitation of the unidirectional rotation for more complicated oblate symmet- 
ric top molecules, like benzene, and compare this analysis with the full quantum-mechanical 
treatment of the problem. 

The outline of the paper is as following. In Sec. II we study classically the rotational 
motion of a linear molecule kicked by a pair of short nonresonant laser pulses that are 
linearly polarized in arbitrary directions. In Sec. Ill, we describe the response of a thermal 
ensemble of linear molecules to such an excitation by using the Monte Carlo method. Angular 
distribution function and its time-averaged behavior are studied in Sec. IV, and the results 
are compared to the full quantum mechanical treatment. In Sec. V, we analyze (both 
classically and quantum mechanically) the problem of inducing unidirectional rotation of 
benzene molecules. Sec. VI summarizes and concludes our study. 



II. CLASSICAL DYNAMICS OF A KICKED LINEAR MOLECULE 

Here we consider classical dynamics of a model diatomic (or, more generally, linear) 
polarizable molecule that interacts with an ultrashort laser pulse. The molecule is treated 
as a linear rigid rotor. We assume that the central frequency of the pulse is far detuned from 
any molecular transition. The potential energy of the pulse interaction with the induced 
molecular polarization is given by j22|: 

V{e, y., t) = -\s\t) (A« cos^ /3 + «x) , (1) 
4 



where Aa = a\\ — a± is the difference between the polarizabihty along the molecular axis 
and the one perpendicular to it, £(t) is the envelope of the electric field of the linearly 
polarized laser pulse, and (3 = I3{6,(f) is the angle between the molecular axis and the 
direction of polarization of the pulse. Here 6, ip are the polar and the azimuthal angles 
characterizing the direction of the molecular axis, respectively. We assume the pulse duration 
is very short compared to the rotational period, so that the pulse can be described in the 
impulsive (5-kick) approximation. We define the dimensionless interaction strength P, which 
characterizes the pulse, as 

A POO 

Let us remind again the pulse sequence producing unidirectionally rotating molecules 



17H19|: (1) at time t=0 molecules are kicked by a pulse linearly polarized along the z-axis; 
(2) at the time of the maximal molecular alignment (t=tah which is the moment when the 
alignment factor (cos^6') reaches the maximum value), a second pulse is applied, polarized 
at 45 degrees to the z-axis (we choose the polarization vector to be in the x-z plane). We 
represent the molecule by a unit vector pointing from the center of mass towards one of the 
atoms (the center of mass is defined as the coordinates' origin). The tip of the vector defines 
a point on the surface of the unit sphere. The orientation of the unit vector is given by 

r = {x,y, z) = {sinO cosip,sm9 sm{p,cos9) . (3) 

The rotational velocity of the unit vector r can be described by two orthogonal angular 
velocities vq = 9 and = {psm9 (the dot denotes a derivative with respect to time), or, in 
vector notation, 

V = veee + v^e^ , (4) 

where the spherical unit vectors are 

60 = {cos 9 cos ip, cos 9 simp, — sin 9) ; 

= (— sinyj, cos(y9, 0) . (5) 

This rotational velocity is referred to as simply "velocity" in the following. The rotational 
kinetic energy in this notation is given by 

T = ll [vl + vl) , (6) 
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where / is the moment of inertia of the molecule. 

As the first step, we find the change in the molecular velocity due to interaction with a 
single short laser pulse. For a 2;-polarized pulse, angle {3 is equal to the polar angle 6*, and 
we can use the potential energy from Eq. ([T]) to find the torque T applied to the molecule 
in spherical coordinates: 

Newton's second law gives Vq = Ivg = 19 and = Iv^ = /J^(0sin6'). We assume that the 
molecule does not change its orientation during the short pulse (impulsive approximation). 
By integrating the Newton's equations of motion over the short duration of the pulse, we 
obtain the change of the velocity due to the laser pulse action: 

At;^ = -yPsin2^o ; Af <p = (8) 

Integrating again, and assuming that the molecule was initially at rest, we get 

eit) = eo-jPtsm2eo ; ifi = ifio. (9) 

Here and (fo describe the initial orientation of the molecule. Every molecule starts to 
rotate with a constant velocity (depending on Oq and ipo ), and r is tracing a circle passing 
through the north and the south poles of the unit sphere. 

As the next step, we consider the action of a pulse linearly polarized along some arbitrary 
unit vector p, and determine the vector of the resulting velocity change Av for a molecule 
oriented along some direction Tq. The norm |Av| can be found similarly to the previous 
case of a 2;-polarized pulse, and it is equal to ||Psin2/3o|, where /3o is the angle between 
the polarization direction of the pulse p and the initial orientation direction of the molecule 
To). The direction of Av can be defined from the vector diagram in Fig. |2l Notice that 
Av is always perpendicular to fq. Also, Av is directed parallel or antiparallel to the vector 
component of p perpendicular to fq, which is equal to p — cos/3oro- As a result, we arrive 
at: 

2hP 

Av = — p cos f3o (p - cos /3oro) . (10) 

We can easily check that ( |T0|) gives the correct result for 2;-polarized pulse, for which p = 
(0,0, 1), by transforming (jH]) to cartesian coordinates using Eq. (jl]). 
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FIG. 2: A linear molecule, represented by ro, is kicked by a laser pulse linearly polarized along an 
arbitrary direction p. As a result, it receives a velocity A v. 



Finally, we provide a formula for the time dependence of the position r and velocity v of 
the point representing molecular orientation on the unit sphere during the field-free rotation, 
given the initial orientation tq and initial velocity vq. As v is always perpendicular to r, 
the tip of the r vector traces a circle on the sphere. The orientation of the molecule at time 
t is 



III. ENSEMBLE AVERAGING AND MONTE CARLO SIMULATION 

We are interested in the response of a thermal molecular ensemble to the double-pulse 
excitation, and in the calculation of ensemble averaged values of the alignment factor and 
other observable quantities. One way to calculate the time dependence of an observable 
averaged over the ensemble is to find its value at time t using (|9]), or ( ITT]) , and to integrate 
it over the initial angular distribution function For example, at zero temperature, the 
molecules are initially at rest and isotropically distributed in all directions. Therefore, for 
any observable f{9, ip), we have to substitute the time-dependent expressions (jH]), ( ITTl) for 9 




(11) 



where vq = |vo|. Taking the derivative, we find the molecular velocity: 




(12) 
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and to obtain f{9o,fo,t), and to integrate over the isotropic distribution of and ipo- 

2tt it 



(13) 





Here the brackets denote averaging over the ensemble. For a nonzero temperature, an averag- 
ing over the initial velocity distribution should be performed as well. In a thermal ensemble 
at temperature T, the velocities vg^ and v^^ are distributed according to the Boltzmann 
distribution: 

Prob(w;_0'<,o)c^^e,o^<,o 
^ : exp 



1 



exp 



"2^ 



dv' 



(14) 



where the dimensionless width of the thermal distribution cxth is given by a^^ = IksT /h^ = 
{2hBc/kBT)~^. Here B is the rotational constant of the molecule, and v' is the dimensionless 
velocity Iv/h. Instead of f|T3|) . we have 

2-K TV OO 

ifir)) = ^ j dipo j dOosmOo j dv'gQ (15) 



dv' Prob(we 0. vLo)fi^o, <^o, Vq q, v' t) . 



These integrals have fast oscillating integrands and are difficult (although possible) to cal- 
culate numerically. 

A somewhat simpler way of finding the time dependence of observables averaged over 
the molecular ensemble is by means of a Monte Carlo simulation. We start with molecules 
isotropically oriented in space, with rotational velocities distributed according to Eq. f|T^ . 
These initial values are obtained using a numerical random number generator. In order to 
randomly distribute points on the surface of the sphere, one uses the probability function 



Prob(^, ^)d9d^ = PToh{9)d9 x PToh{^)d^ 
= -sm9 dO X —dip . 

2 ZTT 



(16) 



The realization of Prob(y9) is straightforward if one has a uniform random number genera- 
tor. If denotes a random number from a uniform distribution in the interval [0, 1], then 



8 




rev 



FIG. 3: (Color online) The ensemble-averaged azimuthal factor (cos^ if) as a function of time after 
the second pulse exciting unidirectional molecular rotation. It is calculated by the classical Monte 
Carlo method with 10^ molecules (dashed black), and by the exact quantum (finite-difference time- 
domain) simulation (solid red). The polarization vectors are pi = (0,0,1) and p2 = (1,0, l)/\/2, 
and P = 5 for both pulses. The calculation was done for nitrogen molecules at the rotational 
temperature of 50 K. The time is measured in units of the quantum revival time Tj-ev = 2TTl/h. 



^0 



rand 



27rw^n- On the other hand, in order to realize Prob(0) one sliould use the transfor- 



mation method [23, section 7.2]. Starting with the equality of probabilities dw = ^ \ sin 9d9\, 
we obtain ^g'^'^'^ = 2 arcsin ^/w^. 

In order to realize the velocity distribution ([H]), one uses the normal random number 



rand 



generator. Using the transformation method again, we obtain v^^q 



O^th^n 



and v^^J^"^ 



(Tth^norm, where Wnorm denotes a random number from a normal distribution with zero mean 
and standard deviation of unity. 

Now we have all the ingredients needed to simulate the scheme for exciting unidirec- 
tionally rotating molecules, as explained at the beginning of Sec. II. Randomly oriented 
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molecules are "kicked" by a pulse polarized along the z-axis, and get velocity additions ac- 
cording to Eq. ffTOj) . which should be added to their initial thermal velocities. The molecules 
then evolve freely according to Eqs. (ITT]) and f ll2p . At each time step, we calculate the 
ensemble averaged alignment factor (cos^O). At the time of the maximal alignment, the 
45-degrees rotated pulse is applied, velocities are changed again according to Eq. ( TTUj) . and 
another free propagation follows. 

As the first application of this approach, we investigate the azimuthal factor (cos^ (p) 
that characterizes confinement of the molecular angular distribution to the x-z plane. This 
factor takes the unity value when all the molecules are confined to this plane, and it is 
equal to 0.5 for the isotropic ensemble. In Figure [31 we plot the classically calculated 
azimuthal factor (cos^ ip) as a function of time following the second pulse, along with the 



exact quantum result (for details of the quantum calculation see 19|). The calculation here 
(and also at Figs. [Hand |5] below) is done for nitrogen molecules with the rotational constant 
of -B = /i/(87r^/c) = 2.00 cm~^. The rotational temperature is 50 K, or (Xth = 2.94. 

It is clearly seen that the classical calculation (Monte Carlo realization using 10^ 
molecules) perfectly reproduces the first oscillation of the azimuthal factor. Obviously, later 
fractional and full revivals are quantum features, which are not reproduced by the classical 
treatment. It is also evident that the time averaged azimuthal factor (that is around 0.52 
here) coincides for the classical and the quantum calculation. This value is higher than the 
isotropic value of 0.5, hence the molecular distribution is squeezed towards the plane defined 
by the polarization directions of the first and the second pulses. We see that a relatively 
simple classical model reproduces the important features of this problem. 

IV. ANGULAR DISTRIBUTION FUNCTION 

In this Section we estimate classically the angular distribution function of an ensemble 
of linear molecules excited by two pulses, and compare it with the results of full quantum- 



mechanical simulation [19|. To achieve this goal, we use the kernel estimator method [24 1 
to reconstruct the distribution function from a set of >> 1 points on the unit sphere, 
which are participating in the Monte Carlo simulation. According to this method, each one 
of the A^ points is surrounded by a narrow symmetric Gaussian function on the surface of 
the sphere. The sum of all Gaussians represents the angular distribution function. Its time 
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dependence is determined numerically by propagating independently rotating molecules, 
as explained in the previous sections. 

A simple way to build a symmetric Gaussian function around a point ri(t) on a sphere 
[i denotes the molecule number in the Monte Carlo realization) is by referring to the angle 
a between the unit vector to this point, and the unit vector to a general point r in its 
neighborhood: cos a = r ■ ri{t). The corresponding (unnormalized) distribution function is 
exp (— Q;^/2cr^), where a is the distribution width (a free parameter). Because only small 
width values are used (cr ^ 1), the quantity can be replaced by the expression 1 — cos a, 
which is equivalent to it for a ^ 1, while for a > 1 the Gaussian distribution function is 
negligible any way. Normalizing the two dimensional surface Gaussian function, we finally 
obtain the following estimate for the time dependent distribution function: 

An important object showing the symmetry and anisotropy of our system is the time- 
averaged angular distribution function. A straightforward, but somewhat cumbersome way 
to generate it is by numerical averaging the distribution function f lT7|) over a long time 
period. A much simpler approach to time averaging is described below. After the end of 
the pulses, each molecule moves with a constant speed on some circle on the surface of the 
unit sphere. Within a sufficiently long time period, each molecule spends the same amount 
of time at each length element of the circle (apart from very rare cases in which a molecule 
stays strictly at rest). Thus, the long-time-averaged angular distribution function for each 
molecule is a "belt" (a circle with a width) on the unit sphere. Its orientation depends on 
the initial position tq and velocity vq on the unit sphere, and its width is a free parameter 
(Tbeit, which, for consistency, is taken equal to the parameter a that was introduced before. 
The total time-averaged distribution function is obtained, as before, by adding the "belts" 
for all the molecules participating in the Monte Carlo simulation. 

The orientation of the "belt" for the i-th molecule is given by the angular momentum 
unit vector cl = ro,i x Vo,i/wo,i, which is perpendicular to the plane defining the circle of 
rotation. Thus, the Gaussian "belt" for a single molecule is located close to the points r 
on the unit sphere, for which the dot product Cl ■ r is nearly zero. Inserting the correct 
normalization factors, we finally obtain for the time-averaged angular distribution function. 
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FIG. 4: (Color online) Time-averaged angular distribution function for an ensemble of nitrogen 
molecules kicked by a laser pulse with P = 10. The curve plots and inset (a) are for molecules 
initially at zero rotational temperature. The dashed green curve presents the exact quantum result, 
and the dash-dotted red line displays the analytical classical distribution according to Eq. (fT9|) . 
The blue solid curve and inset (a) are calculated by the Monte Carlo method according to Eq. 
(jlSp . with cJbeit = 0.1 and N = 10^. Inset (b) shows Monte Carlo results for the initial rotational 
temperature of 50K. 



Pta- 



(18) 



1 ^ 



'belt 



exp 



where the overline denotes a long time average. 

Before considering the double pulse scheme, we first analyze the time-averaged angular 
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distribution function for molecules subject to a single pulse polarized along the z-axis. The 
result for zero temperature (all molecules are initially at rest) is shown in Fig. H] (a) and 
duplicated at the plot by the solid line. In Fig. |l](b), the distribution function is plotted for 
the same system, but at finite temperature. As expected, the distribution is "cigar" -shaped, 
and azimuthally symmetric. At zero temperature, an analytical result can be obtained for the 
time- averaged distribution function in the classical limit. In this case, all the trajectories of 
the kicked molecules are circles passing through the poles of the unit sphere. These circles 
are uniformly distributed in the azimuthal direction, and the distribution function is ip- 
independent. Therefore, at any given value of 6 we find the same total number of molecules. 
As a result, the time- averaged distribution function multiplied by the geometrical factor 
sin 6' should be a constant. It is seen from here, that the distribution function averaged over 
long time should be proportional to 1/ sin^, or, after normalization: 



This equation was verified also by evaluating numerically Eq. (9) of reference 25|, including 
summation over all branches. The function of Eq. (fT9|) is plotted in Fig. lUby a dash-dotted 
line. Our numerical Monte Carlo result in Fig. |4] (solid line) agrees well with the analytical 
result, but it is not singular at the north and the south poles, because a finite width abeit 
was given to every molecule trajectory, and the singularities are smeared out. 

Finally, we compare the classical result to the quantum angular distribution function 
averaged over a single revival period. The quantum result is given in Fig. H] by a dashed line. 
Two features should be stressed here. First, the singular classical behavior disappears in the 
quantum description. Second, there is a fine oscillatory structure in the quantum distribu- 
tion, which is absent in the classical distribution. Besides these features, the distributions 



are similar. We also compare the values of (a;^), {y^) and (-2^), calculated both classically and 
quantum mechanically at the conditions shown in Fig. H] (b). Here the brackets denote the 
ensemble averaging and the overline denotes long-time averaging. The quantities x, y, z are 



the cartesian components of the unit vector r. In both cases, the same values of (2;^) = 0.42 



and (x^) = (y^) = 0.29 are obtained. This should be compared to the "isotropic" values of 



(a;2) = = (^2) = 1/3. 

For the double pulse scheme, leading to the unidirectional rotation of the molecules, we 
present the classical time-averaged angular distribution in Fig. |5] (a), which is very similar, 
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FIG. 5: Time-averaged angular distribution function for unidirectionally rotating nitrogen 
molecules that were excited by two delayed and cross-polarized pulses. The initial rotational 
temperature is 50 K. The calculation was done using the classical Monte Carlo method, Eq. (jlSp . 
in (a), and using the quantum finite-difference time-domain simulation followed by the averaging 
over a revival period in (b). Both pulses had the kick-strength of P = 5. The following parameters 
were used in (a): fibcit = 0.1, and N = 10^. 

except of fine details, to the corresponding quantum distribution in Fig. |5](b), presented also 
in [l9|. The distribution is squeezed to the xz-plane (the plane defined by the polarization 
directions of the two pulses). We expected (y^) to be smaller than (x^) and {z'^), and we 
indeed obtain numerically (x^) = 0.31, (y^) = 0.30 and {z"^) = 0.39 classically, as well as 
quantum mechanically. These results demonstrate very good agreement between the classical 
and quantum calculations even for relatively small values of the interaction strength P. 
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V. UNIDIRECTIONAL ROTATION OF BENZENE MOLECULES 



Now we address a more difficult problem of the unidirectional rotation of benzene 
molecules. These molecules, considered as rigid bodies, behave as oblate symmetric tops. 
The classical treatment of the free rotation of a symmetric top is analytically simple (al- 
though more complicated than the rotation of a linear rotor), and is well documented in 
literature (see, e. g. [26, §§33,35]). 



A. Rotation of benzene molecules — classical treatment 



There are two moments of inertia characterizing a symmetric top, Ji = I2 and I3, which 
are defined using the body fixed coordinate system xi,X2,Xs, shown in Fig. [6l where the 
molecular symmetry axis is directed along the X3 axis. The figure is taken from 26|, §33]; 
the X2 axis is perpendicular to the plane of the figure and is not shown. For benzene 
molecule, which is an oblate symmetric top, and has the form of a planar ring, the X3 axis 
is perpendicular to the plane of the ring and passes through its center. Hence, for benzene. 



Ji + /2 = 2Ji = h 



m §32]. 



When the molecule is freely rotating, the angular momentum, denoted by L, is constant 
in magnitude and direction. The rotational kinetic energy is also a constant of motion and 
can be written as: 

2/1 2/1 2/3 2/1 2J3 ' ^ ^ 

where Li, L2 and L3 are the components of angular momentum along the body- fixed coor- 
dinate axes Xi, X2 and X3, respectively, and L| = + L"^. The free rotation is such that 
the angular momentum L, the angular velocity O and the molecular axis X3 all lie in the 
same plane. The molecular axis X3 (as well as the vector ft) performs precession around the 
angular momentum vector L with an angular velocity Opr, which is shown in Fig. |6l Its 
magnitude is: 

^r>r = T ■ (21) 
-'1 

The molecular axis 0:3 moves on the surface of a cone and forms an angle 6'pr with the cone 
axis (which is along the vector L). 

In addition to the precession motion, the top rotates around its X3 axis, but this motion 
is not important in our problem. The reason is that this last rotation does not change the 
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FIG. 6: The free rotation of a symmetric top with notation used in the text. 

molecular orientation (defined by the molecular axis), and that the laser pulse does not 
influence this motion, as will be explained later on. 

Wc describe the orientation of the molecular axis by a unit vector r in the direction of 
x^. Its value at t = is Tq. The velocity of this vector is given by 

V = ripr X r . (22) 

We now define two more unit vectors. One is cl, a unit vector in the direction of L. Another 
is in the plane of the circle traced by the tip of r: 

ro-cosgp,eL 
^ sin^p, ■ (2^) 

Using the notation defined, it is easy to see that the time dependence of r is given by: 

r(t) = cos ^prBL + sin [ roll cos Qpri + — sin f2prt | (24) 

For the velocity we obtain: 

v(t) = fipr sin ( — i"o|i sin Vt^j-t H — - cos Vt^^t ) (25) 

It is important to stress, that the tip of the unit vector r does not necessarily trace a circle 
of a unit radius on the surface of the unit sphere, as for the linear molecule. It can trace 
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any circle, which can even shrink to a point, in the case of f2i = = and Q3 ^ 0. A unit 
circle is traced only if 6'pr = 7r/2, or 1^3 = and Qi = L/Ii. 

Although we are dealing here with a more complicated molecule, the interaction with an 
ultrashort laser pulse is described by the same formula as in Eq. ([T]). The reason is that the 
polarizability tensor of the symmetric top has the same symmetry as the one of the linear 
molecule. The interaction depends, as before, on the electric field envelope squared, on Aa, 
which is negative in the case of benzene, and on the angle between the pulse polarization 
direction p and the molecular axis orientation fq. Notice that a linearly polarized laser pulse 
does not produce a torque parallel to the molecular axis. 

The torque exerted by the pulse is, therefore: 

r = -^£^{t)Aa sin 2^90 epxro , (26) 

where Cpxro is a unit vector in the direction of the cross product p x tq. This equation 
shows that the laser pulse "kicks" the symmetry axis of the benzene molecule towards the 
plane, which is perpendicular to the polarization vector p. Using Newton's second law 
and integrating over the short time of the pulse, we obtain for the change in the angular 
momentum: 

AL, = -^Psin2/3o(epxro), , (27) 

and two similar equations for the y- and the z-components. Angle Po can be found using 
p ■ To = cos /3o- Now, after the new angular momentum vector is known, and assuming 
the molecule did not move during the interaction with the pulse, the molecular rotation is 
completely defined by Eq. 1^^. Naturally, new values of 6'pr, cl, ro||) vq and Qpr should be 
found using the new value of L. Specifically, 6'pr is found using the relation cos^^pr = cl ■ tq, 
roll is then defined using Eq. (12^ . Qpj. - by Eq. (jH]) and vq - using Eq. fl22]) . 

The Monte Carlo simulation for benzene molecules has several differences compared with 
the one used for linear molecules. Here, the kinetic energy is more easily separable when 
expressed using the angular momentum components along the body-fixed coordinate system, 
as in Eq. (|20|). 

We start with random generation of the angular momentum components Ly and L3, 
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according to the thermal distribution: 



Prob rfLji X 



Prob (L'r^) dL'^ 



(28) 




1 



exp 




X 



1 



th,3 



exp 




dL'^ 



where we introduce a^-^^ = IikBT/h? and a^j^g = I^kBT/h'^, similar to Eq. f|T^ . and where 
L' represents the dimensionless angular momentum L/h. Employing the transformation 
method, and using the notation introduced after Eq. ( |T6i) . we obtain explicitly for the random 
components: (Lj,)''^"'^ = crth,iA/ln [1/(1 — Wun)] and {L'^Y^^'^ = crth.su^norm- Knowing these 



components, we deduce the magnitude of the angular momentum L = \J L'^\ + -^1, but its 
direction is yet to be defined. We obtain also the precession velocity using Eq. fl2T|) . and O^t- 
from the relation cos^pr = L-^/L. 

Next, we determine the direction of the angular momentum in the laboratory coordinate 
system, which is given by the vector cl. We do this by randomly generating two spherical 
angles 6'l and ip-^,. These two angles are distributed the same as 6 and ip in Eq. f ITB]) . because 
the tip of the vector cl is distributed uniformly on the surface of the unit sphere. 

Finally, by randomly generating an angle y?, uniformly distributed between and 2tx, we 
obtain the initial orientation of the molecular axis on the cone around the vector L. Now, 
the coordinates of ro are completely defined in the frame in which the z-axis is fixed along 
the vector L. These coordinates are given by (sin 6'pr cos yj, sin 6'pr sin cos6'pr). Using two 
rotational matrices, the coordinates of fq in the laboratory frame are then found. 

The above Monte Carlo procedure was used to generate data that are discussed later in 



B. Quantum mechanical treatment 

Here we outline the quantum mechanical description of unidirectional rotation of benzene 
molecules by two delayed ultrashort laser pulses with crossed linear polarizations. Details 




Sec. C. 



will be presented in 27|. The molecular Hamiltonian, Hq, in the field-free conditions is 
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obtained simply by replacing the components, Li {i = 1,2,3), of the classical angular mo- 
mentum by the corresponding quantum mechanical operators, Jj, in Eq. (|20|) . The rotational 
eigenstates and eigenenergies are given by the time-independent Schrodinger equation, 

Ho\r) = Er\r), (29) 

where r stands for an index to identify the eigenstates, explicitly expressed by a set 
{J, K, M}, with J being the total angular momentum, and K and M being its projec- 
tion onto the molecular symmetry axis and the space-fixed axis, respectively. At finite 
temperature, an initial ensemble is a mixture of molecules in many different eigenstates, but 
we first take a single rotational level as an initial state, |rj). The nonadiabatic interaction 
with the nonresonant ultrafast laser field, given by Eq. (Q, converts the stationary state to 
a rotational wave packet, |^I/(t)), which is expanded as 

l^n(t)) = U{t,0)\ri) =J2CrAt)eM-^^rt)\r) . (30) 

r 

Here, U{t2,ti) is the time-evolution operator from time ti to ^2 and ojr = Er/h. The expan- 
sion coefficients, Cn^rit), vary during the interaction with the laser pulses, and subsequently 
become constant after the laser field vanishes. They can be derived by solving the time- 
dependent Schrodinger equation (TDSE) with the initial condition, |\E'()f: — )■ — oo)) = \ri), 

^n(t)) , (31) 

where V{t) is given in Eq. ([1]). Substituting Eq. (130|) . this TDSE is recast into the following 
coupled differential equations for the expansion coefficients: 

ihjCrUt) = Y.^AV{t)\r') exp(-zAa;,,,,t)C,,y(t) , (32) 

r' 

with Awr'.r = ^r' —(jJr- Thcsc couplcd equations can be solved numerically to determine the 
complex expansion coefficients, once the matrix elements for V are specified. 

If another laser pulse is irradiated onto the molecule at t = r, the wave packet created 
by the first pulse, represented in Eq. (l30l) . is further modified by the interaction with the 
second pulse. The resultant wave packet at t after the second pulse is expanded as 

\^>rXt)) = U'{t,T)\m,Xr)) = Y.B,,^,{r)eM-^^rt)\r) , (33) 



^0 + V{t) 
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where the prime on U indicates the properties of the second pulse. The quantity -Br^.r is the 
transition amphtude from the initial \ri) state to |r) by the interaction with the two pulses, 
and is represented by 



Br^^riT) = Cr^yCl,^^ exp(-Z AcJ^/ .^r) 



(34) 



To evaluate the transition amplitude given in Eq. flMj) . it is convenient to adopt the 
space-fixed axis system defined differently from the previous one. Here, we set the z axis 
along the laser propagation direction and the polarization of the first pulse parallel to the x 



axis. Then, the interaction given in Eq. ([T]) is recast as 



or, re-expressed as 



V{9,^,t) = -S^{t) {Aa[l + cos{2<^)]cos'^e 
8 

— Aa cos(2(y9) — (a;|| + q;_|_)} , 



(35) 



+ 



D 



(2)* 
-2,0 



+ D 



(2)* 
2,0 



(36) 



where Dp^q is the rotational matrix 



29| . With this expression, the evaluation of the matrix 
elements for a symmetric-top basis set is straightforward by applying the spherical tensor 
algebra, and the following selection rules are derived: 
AK = 0, 

AM = (then AJ = 0, ±2 for KM = and AJ = 0, ±1, ±2 for KM 0), or 
AM = ±2 (then AJ = 0, ±2 ioi K = and AJ = 0, ±1, ±2 for K ^ 0). 
We set the polarization of the second pulse tilted against that of the first one by angle 
A(f. When the new axis system is defined so as the x' axis is set parallel to the polarization 
of the second pulse, the new and old angular coordinates are related as: (p' = ip — Aip. 
Consequently, the symmetric-top wave function transforms as 



I J, K, M)^\J, K, M) exp(iMAv;) . 



(37) 

By taking the transformation into account, B^^^r in Eq- dM]) is rewritten, if the second pulse 
is a replica of the first one, as 



Bn,r{r) = ^ Cr„r'C^',r exp [{-iAur',rr - M'Alf)] 



(38) 
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with r' standing for a set of {J\K\M'}. The transition amphtude for the two pulses can 
be evaluated for an arbitrary r at once, if the set of the coefficients, Cr,,r', are derived by 
solving the coupled equations in (152]) . 

The expectation value of a given observable A for a certain initial condition with finite 
temperature is derived as the following ensemble average, 

(i)(r,t) = Y.WrA^rMM'^rM (39) 

= J2 Wr^B:^y{T)Br^^r{r){r'\A\r)exp{iAur',rt) , 

ri,r,r' 

where Wr^ is the Boltzmann factor for |rj) with appropriate degeneracy factor due to nuclear- 



spin statistics [30|. For the discussion of the molecular alignment and orientation, we take 
cos^ J^, and (corresponding to the classical Ly), as A. 

C. Classical vs quantum treatment 

In this Section we analyze the process of exciting unidirectional rotation of benzene, 
following the setup of [21I. The scheme operation is practically the same as described after 
Eq. ([2]), except of several details specific for the benzene molecules. Shortly after the first 
ultrashort pulse, polarized along pi = (0,0,1), the molecules experience a transient anti- 
alignment, as opposed to alignment in the case of linear molecules. This happens because Aa 
is negative for benzene, as opposed to linear molecules, such as nitrogen, and the direction 
of the angular momentum transferred to the molecule by the pulse is reversed. Therefore, 
shortly after the first pulse the molecular axis angular distribution is squeezed to the plane 
perpendicular to pi. Quantitatively, after the first pulse, the ensemble average of {cos^ 6) 
decreases below the isotropic value of 1/3. The second delayed and tilted laser pulse brings 
this anisotropic ensemble to the rotation with a preferred sense. In the following, we shall 
consider the time development of the alignment factor (cos^ 6) after the first pulse, and the 
dependence of the average value of the induced angular momentum as the function of the 
delay of the second pulse. 

We show first the classical results. The upper part of Fig. [7] presents the alignment factor 
(cos^ 6') as a function of time passed after the first pulse (polarized along pi = (0,0, 1)), 
for three values of the pulse strength P. Here and below we use the following parame- 
ters for the benzene molecule 3l|]: rotational constants B = /i/(87r^/ic) = 0.190 cm~^, 
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2 3 4 

time (units of I^/h) 



FIG. 7: (Color online) Classical treatment. The upper panel shows thermally averaged alignment 
factor (cos^ 9) for benzene molecules kicked by a single pulse, plotted as a function of time after the 
pulse. The pulse has a strength of Pi = —1, —3, —10 (black, blue and red lines online, respectively) 
and polarization direction pi = (0,0,1). The lower panel displays the (normalized) value of the 



oriented angular momentum {Ly) / {!?■) induced by a second ultrashort pulse, shown as a function 
of the delay of this pulse. The polarization direction of the second pulse is p2 = (—1,0, l)/\/2 (at 
45 degrees to pi), and the value of Pi is the same as that of Pi, with the same color code used 
for the curves. The initial rotational temperature is 0.9 K. The calculations were done using the 
Monte Carlo method with 10^ realizations. 



C = /^/(Svr^/gc) = B/2 and the polarizabilities ay = 6.67 A , a± = 12.4 A . All plots 
are for rotational temperature of 0.9 K (as in experiment 2l|), which corresponds to the 
dimensionless parameters crth,i = 1-29 and crth,3 = V^(Tth,i = 1-82. 

The plots for (cos^ 6) show that the anti-alignment becomes more pronounced and occurs 
at earlier times as the pulse intensity increases, as is naturally expected. Similar to the 
case of linear diatomic molecules, we anticipate that the second pulse should be fired at 
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FIG. 8: (Color online) The same quantities as in Fig. [7] are plotted here. They were calculated 
using the quantum mechanical treatment. Half and full rotational revivals are seen, which are not 
present in Fig. [71 

the moment of the maximal squeezing of the molecular angular distribution, i. e. when the 
alignment factor takes the minimal value. In this case, the maximal average torque is applied 
to the molecules, facilitating the most emphasized unidirectional rotation. This is indeed 
seen from the lower part of Fig. [71 Here the quantity {Ly) / ^/{LP) (calculated after the second 
pulse) is plotted as a function of the time delay between the first pulse (with polarization 
direction pi) and the second pulse with polarization direction p2 = (— 1, 0, l)/v^, where 
the angle between pi and p2 is 45 degrees. The maximal modulus of the transferred angular 
momentum happens indeed when the second pulse is applied at the moment of the maximal 
anti-alignment. We also notice that the optimal delay becomes shorter with the strength of 
the pulses. 

Figure [8] presents the same quantities calculated fully quantum mechanically, following 
the procedure outlined in section V-B. On the long-time scale, it demonstrates quantum 
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FIG. 9: (Color online) Comparison of the classical results of Fig. [7] (solid lines), and quantum 
mechanical results of Fig. [8] (dashed lines) at short times. The values of P are —1, —3, —10 for 
black, blue and red lines, respectively. 

rotational revivals, such as half and full revivals, that are absent in Fig. [7] for obvious 
reasons. For a detailed comparison with the classical results at shorter times (corresponding 
to the conditions of the experiment 2l|]), Fig. |9] displays the classical and the quantum 
graphs overlapped on the same plot. The results seem to be practically identical at short 
times, and the agreement is improving with the increase of |P|, as more rotational states are 
involved in the dynamics of the wavepacket. Parameters of the experiment |2l| correspond 
to the curves for P = —3 in Fig. [9], with the delay time chosen near the minima of the 
curves. As follows from Fig. [9l the rotational dynamics in this range is adequately described 
by the classical treatment. 
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VI. CONCLUSIONS 



We studied the problem of exciting unidirectional molecular rotation with the help of two 
delayed cross-polarized laser pulses using classical mechanics, and compared this approach 
with the results of fully quantum mechanical treatment. The problem was analyzed both for 
linear molecules and for symmetric tops, like benzene. Although the comprehensive descrip- 
tion of the fine details of rotational dynamics undeniably requires the full scale quantum 
formulation, we found a very good agreement between the classical and quantum approaches 
in a wide range of experimentally relevant parameters. The classical model identifies the 
origin of the mechanism, defines the optimal delay between the pulses, and the best an- 
gle between them. Moreover, it accurately predicts the magnitude of the induced oriented 
angular momentum, and correctly reproduces the behavior of the optimal delay with pulse 
strength, as it was observed in 2l| . At the same time, a fully quantum treatment is required 
to describe quantitatively the operation of the double-pulse scheme for relatively weak laser 
pulses, when only low rotational states are involved. In this case, considering quantum inter- 



ference of excitation pathways by the first and the second pulses [2l|] provides an instructive 
complementary view on the origin of the phenomenon. Moreover, operation of the scheme 
at long delays between the pulses (comparable with the period of rotational revivals) needs 
quantum description as well. At the same time, quantum and classical results agree well 
again if the long-time-averaged observable quantities are considered. We illustrated this by 
comparing classical and quantum time-averaged molecular angular distributions that show 
anisotropy and confinement to the plane spanned by two polarization vectors of the laser 
pulses. 

We believe that this comparative study provides the reader a comprehensive view on the 
physics behind the double pulse scheme for ultrafast preparation of molecular ensembles 
with oriented rotational angular momentum and controlled sense of rotation. 
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